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In this paper, we describe several new techniques for use in the 
design of switched communications networks. These techniques apply 
to the development of traffic routes which realize network traffic 
flows in the context of an existing optimization method that assigns 
these flows. The general ideas involve the careful selection of basic 
variables and the successive reduction of the problem to one of convex 
hull formation in Euclidean n- space and finding Hamiltonian cir- 
cuits for a class of highly structured graphs. We include several 
examples showing how these techniques are applied. 

I. INTRODUCTION 

Recently, R. H. Cardwell 1 proposed a switched communications 
network design algorithm for the future stored program control net- 
work. The networks under consideration are nonhierarchical in struc- 
ture and take advantage of traffic noncoincidence in routing. The basic 
objective of Cardwell's algorithm is to design a minimum cost trunking 
network which, by using an appropriate routing strategy, can carry the 
necessary traffic load and, at the same time, meet the required grade 
of service. 

In this paper, we describe an extremely efficient method for produc- 
ing an appropriate routing strategy. One of our original intentions was 
to develop a mathematical framework into which dynamic routing 
problems, such as those described later, could be placed. Indeed, it 
seems likely that the approach used here may be valuable for exam- 
ining other classes of such routing problems. 

II. BACKGROUND 

In Fig. 1 we show a block diagram of Cardwell's algorithm. Suppose 
we wish to design a network using the algorithm. (See Ref. 1 for a 
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Fig. 1— Block diagram of Cardwell's algorithm. 



more detailed description.) We start by initializing the blocking prob- 
abilities of each link. The routing module selects a set of the most 
economical paths for each pair of nodes and then assigns flow to the 
paths. Routes, which are ordered lists of paths, are then formed so 
that the probability that all paths in any list are busy is small enough 
to meet the required grade of service. Then, by means of a linear 
programming formulation, the routing module determines a network 
flow which minimizes the total cost, considering link costs and traffic 
noncoincidence. In the engineering module, the Erlang loss formula is 
then used to fix the number of trunks required for each link. 2 In the 
update module the ECCS method of Truitt is used to help minimize 
the network cost. 3 The blocking probabilities for all links are then 
updated. The whole process is now iterated until satisfactory conver- 
gence is achieved. 

Figure 2 shows a block diagram of the routing module for the unified 
algorithm. A basic feature of our method is that actual routes are not 
formed until convergence has been obtained in an earlier part of the 
unified algorithm. Only after this occurs does the routing realization 
submodule generate the routes and provide the appropriate routing 
strategy. Refer to the work of Murray and Wong which gives efficient 
heuristic algorithms for solving the linear programming problems in 
this module. 4 The upper bound module is a new addition which helps 
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Fig. 2— Modified routing module. 
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the iterative procedure to converge more rapidly by setting stronger 
upper bounds for the carried loads of the various paths chosen. 

One of the questions concerning this algorithm was the problem of 
synthesizing routes from the assigned path flows (the route realization 
block in Fig. 2). The following is a discussion of an efficient technique 
that accomplishes this. 

III. CYCLIC ROUTING 

Consider the special case of two nodes A and B. Assume there are 
n paths Pk, 1 S k < n, between A and B. The amount of traffic to be 
carried on P* is denoted by Xk, where we have normalized the traffic 
load so that one unit of traffic is attempted between A and B. The 
blocking probability for Pk is denoted by p*. We will call the vector 
x = (xi, ■ • • , x„) the desired traffic vector and p = (pi, • • • , p„) the 
blocking probability vector. (The x's are actually outputs of the linear 
programming module.) 

For a permutation m of {1, 2, • • • , n), by the route R{tt) generated 
by tt, we mean the route in which the path P w (d is tried and, if blocked, 
path P„(2) is tried. If that path is blocked, then path P„m is tried, etc. 

The first question is: what are the traffic flows on the various Pk 
when route R (tt) is used? Let qk = 1 — Pk and assume that tt is the 
identity permutation, i.e., 7r(k) = k for all k. Since one unit of traffic is 
initially attempted on Pi, the first path of R (tt), then q t units of traffic 
are carried on Pi and pi units of traffic are blocked. These pi units are 
now attempted on P2. Thus,pi#2 units get carried andpip2 are blocked. 
Continuing this process, we see in general that on Pk, 
P1P2 • • • Pk-\qk units of traffic are carried and P1P2 • • • Pk-iPk are 
blocked. We condense this information into the flow vector ¥(tt) = 
(Fi(tt-), F 2 (tt), • • • , F„(tt)) = (<7i, piq 2 , PxPiq.u • • • , P1P2 • • • q n ). Note 
in particular that the amount of traffic which is blocked is just 
pip 2 • • • p„, independent of tt. 

The overall plan is to use each route R (tt) a certain fraction acin) of 
the time, as tt ranges over all permutations of (1, 2, • • • , n), so as to 
achieve the desired traffic flow Xk on each Pk. In other words, if 
possible, find a (tt) with 

a{ir) > 0, X oc(tt) = 1 
so that 

X = ]>>(77) F(tt), 

where tt ranges over all permutations of (1, 2, • • • , n}. 

This is exactly the same problem as deciding whether x is in the 
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convex hull of points F(7r) (considered as points in n-dimensional 
Euclidean space E n ) and, if so, finding a representation of x as a 
convex combination of the F M . Note that all the F(7r)'s are extreme 
points of the convex hull. Since 

£ Fa(tt) = 1 - pxpi • • • Pn for any tt, 
k 

then the convex hull is actually (at most) an (n — 1) -dimensional 
polytope. Thus, any point in the convex hull can be represented as a 
convex combination of some choice of n extreme points ¥{ir). 

As an example, we consider in detail the case n = 3. In Table I, we 
list the six possible tt's and the corresponding F(tt)'s. 

We will denote the permutation it which sends i to tt(i') by the 
sequence ir(l)ir(2) • • • tt{jx). This should not be confused with the 
ordinary cycle notation for a permutation tt (which will also be used). 
For example, the permutation of {1, 2, 3, 4, 5, 6} given by tt(1) = 3, 
tt(2) = 5, tt(3) = 6, it{4) = 4, tt(5) = 2, tt(6) = 1 can be written both as 
•n = (136)(25)(4) and it = 3 5 6 4 2 1. 

Figure 3 shows a typical picture when these points are plotted in E 3 . 
All six points he on the plane Fj + F 2 + F 3 = 1 - P1P2P3. We should 
note here that we always assume < p k < 1 for all k, since any path 
with blocking probability one can be removed without affecting the 
traffic flow, and any path which carries any traffic at all has positive 
blocking probability (less than one). 

In general, we would like to be able to decide if the desired traffic 
vector x lies in the convex hull of F(tt) and, if so, how to represent it 
as a convex combination of F(tt). A natural choice to consider is a 
cyclic set of routes. For example, suppose we consider the three routes: 

7T, = 1 2 3, 

7T2 = 2 3 1, 

7T3 = 3 1 2. 

Let us determine whether x is in the convex hull of these three points. 
Of course, a necessary condition is £,• x, ■= 1 — pt p2pa- In any case, 

Table I — Flow vectors for the 
Case n = 3 
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F : +F 2 + F 3 - 1-p,p 2 p 3 




Fig. 3— Geometrical representation of the flow vectors for the Case n = 3. 

since the convex hull is 2-dimensional, any point in it is a convex 
combination of some set of three extreme points. The cyclic sets seem 
reasonable choices since they apparently span rather large portions of 
the convex hull, although certainly not all of it. For example, in Fig. 3 
we have shaded the convex hull of F(tti), FM, and F(ir 3 ). This is 
much larger than, say, the triangle spanned by F(123), F(132) and 
F(231). 

Therefore, we are looking for coefficients a, such that 

£ <*Ffo) = x, (1) 

1=1 

with 

a, > 0, £ a, = 1. 
i 

By eq. (1), the a, must satisfy 

I euFkte) =x k ,k = 1, 2, 3. 

i-l 

Expanding these equations using Table I, we obtain 
ot\q\ + a2<7i/?2/>3 + otaqipa = Xi, 
ot\P\qi + OL2q2 + cop^zpa = *2, 

a\P\p-iqA + cupzq* + a.3^3 = X3. (2) 
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The determinant A of the system eq. (2) is given by 

<7l QlP2P3 QlP3 

A = p\Q2 Q2 P\q^qz 
PxPiqz P2?3 qa 

1 P2P3 pa 
= qiq2qa Pi 1 P\Pa 

PlP2 P2 1 

= 919293(1 - P1P2P3) 2 . 

Solving for the a*, we have 

Xi q\P2P3 qiP3 



«, = - 



x 2 92 Piq2qa 
X3 P2qs 93 



Letting 



= [(xi/qi) - (p3X 3 /q3)V(l ~ P1P2P3), 
(X2 = [(X2/92) - (piXi/qi)]/(l - P1P2P3), 
as = [(xs/qs) - (P2x 2 /q2)]/(l ~ P1P2P3). 



Xk . PkXk 

Ok = — , Ok , 

qk qk 



we see that the a* are > if o\ — S3 > 0, 02 — 5i ^ 0, 03 — 82 S 0. 

For general n, a similar calculation shows that the corresponding 
system of n equations has determinant A given by 

A = q t q 2 ■ • • q n (l - P1P2 • • • p n ) n ~ l 

and coefficient values 

a*+i = [(** + i/9* + 1) ~ (Pkx k /qk)]/(l ~ P1P2 ••• Pn) 

for the cyclic set of routes 

1234 n 

234 .... n\ 
3 4 .... n 12 

n 1 2 ..'•• n- 1 
where addition of indices is modulo /i, i.e., 

a= (01- 8 n )/(l - p\ ••• Pn). 

Consequently, we succeed with this cyclic choice of routes if all the 
ak's are at least 0, i.e., 
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02 > Si, <J3 > 82, • • • , 0„ > S n -1, 0\ > S„. (3) 

Note that £,- «< — 1 follows at once from 

£ jfi f - 1 - p x • - • p» 

i 

and, in particular, note that the labeling of the Pk is arbitrary. Any 
arrangement of the a's and S's satisfying eq. (2) will give us a cyclic set 
of routes which works, i.e., a set of routes which contains x in the 
convex hull. 

In order to find these efficiently, we can do the following: From the 
given Xk and pk form 

qk = 1 — Pk, 

x k 
Ok = — , 
Qk 

« PkXk 

Ok = . 

Qk 

We are just searching for a cyclic permutation (j\, jz, • • • , jn) of 
{1, 2, • • • , h) such that 

q, 2 > 8,,, o/ 3 > 5y 2 , • • • , a,, > 5/ n . 

To find this, form the directed graph G which has as its vertex set the 
set of paths Pk and an edge from each P, to Pj for which q, > 8». If in 
G we find a Hamiltonian circuit (i.e., a circuit passing through each 
vertex exactly once), say, Pj { Pj 2 • • • Pj n {Pj), then by the definition of 
the edges of G, we must have 

Oh ^ *>,. a /> s $/ 2 » ' ' * • °A - 8 Jn> 

which is precisely what we want. 

Thus, we have shown that x can be realized from a cyclic set of 
routes if and only if G has a Hamiltonian circuit. Of course, the 
problem of finding a Hamiltonian circuit in an arbitrary graph is 
known to be an NP-complete problem (see Ref. 5 for an exposition of 
this term) and, therefore, almost certainly computationally intractable 
as the graph becomes large. Fortunately, however, the graphs G are 
far from arbitrary and, in fact, we can provide an algorithm for finding 
Hamiltonian circuits in them which runs in time 0(n log n). 

First, we may assume without loss of generality (by a suitable 
relabeling) that 

o\ > a 2 > • • • > o n . 

Note that a necessary condition for the existence of a Hamiltonian 
circuit in G is: 
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For all k, 2 < k < n, 

|{i:a*>5,}| >n-k + 2, (4) 

where \X\ denotes the cardinality of the set X. To see this, note that 
if G has a Hamiltonian circuit, then for each k, there is at least one 
edge from a vertex in {P k , Pk+i, ••• , Pn) to one in {Pi, P 2 , - • • , 
P A _i}.Thus, 

Ok ^ or > S t " 

for some t', t" with t' > k > t". Therefore, 

{i:a*>fc}D{i:i>*} U {*"}, 

which implies 

| {i.ok >8i}\>n-k + 2. 

In fact, eq. (4) is also a sufficient condition for G to be Hamiltonian. 
This can be seen from the following proof (by induction on n). 

Suppose n = 2 and eq. (4) holds. Then clearly a 2 > 5i and G is 
Hamiltonian. Next, assume that eq. (4) is sufficient for all such graphs 
with n — 1 vertices. Suppose G has n vertices and satisfies eq. (4). Let 
G' be the induced subgraph on {Pi, P 2 , • • • , P*-i} (where we have 
assumed as usual that ai > a 2 ^ • • • > a„). It is easy to see that G' also 
satisfies eq. (4). By induction, G' has a Hamiltonian circuit, say, 
Pj Y Pj 2 • • • P/ n _,- Since G satisfies eq. (4), 



But also 



Thus, 



a n S 8j. for some i, 1 < i < n — 1. 

a k ^ On ^ 8 n for all k, 1 < k < n - 1. 
P . . . P P P . . . P 



is a Hamiltonian circuit in G and the induction step is completed. This 
proves that eq. (4) is, in fact, a necessary and sufficient condition for 
G to have a Hamiltonian circuit. 
We summarize the preceding discussion in the following. 

Algorithm for cyclic routing 

Inputs: Paths Pi, • • • , P n joining two given points A and B, and the 
corresponding blocking probability vector/? = (pi, • • • , p n ) and desired 
traffic vector x = {xi, • • • , x n ). 

Object: To find a permutation ir of {1, 2, • • • , n) satisfying 

i-l 
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with 

a, > 0, £ a, = 1, 

where F is the flow vector function and in is the cyclic route i, ir(i), 
7r {2) (i), • • • , w (B_1) (i) (i-e., P« is tried first, then P w(l) , etc.). 
Algorithm: 
(i) Calculate 

a* = and 5* = for 1 < & < n. 

1 - p h 1 - Pk 

(Recall that we are always assuming that < pk < 1 for all k .) 

(ii) Relabel the a*, if necessary, so that o\ > 02 > • • • ^ a n . 

(Hi) Set it <- (1), i <r- 2, y «- 5 lf z «- 1. 

(w) If a, < ^, go to (i>«). If a, > % insert i after 2 in the cycle 
representation of it. 

(v) If v > Si, set v <— 5„ 2 <— i. If y < 6«, y and z are unchanged. If 
i < n, set i <— i + 1 and go to (iv). 

(vi) The desired Hamiltonian circuit is P\P„ { i)P„m {i) • • • P„«p-i)(i). 
Define 

Ok — S„-\(k) , 

a* =- , 1 < ^ < n. 

1 - P1P2 • • • Pn 

£ can be realized by using route tth for the fraction a* of the time, 1 < 

(um) x cannot be realized by any set of cyclic routes. 

End. 

Note that except for (ii), in which n log 2 n operations are required 
in the ordering of the n 0,'s, all other steps require at most 0(n) 
operations. Thus, the computational complexity of the algorithm is 
n log2 n + O (n ) in time and O (n ) in space. 

We point out that the desired traffic flow vector x can often be 
realized by more than one set of cyclic routes, i.e., the graph G might 
have more than one Hamiltonian circuit (each of which corresponds 
to a cyclic routing realization). The preceding algorithm will always 
produce one such realization provided any exists at all. 
Two examples 

Example 1: There are five paths between two points A and B. The 
desired traffic vector x and the blocking probability vector p are as 
follows: 

x = (0.185, 0.231, 0.220, 0.072, 0.242) 

p= (0.8,0.7,0.6,0.5,0.3). 
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Thus, 



a = (0.924, 0.770, 0.550, 0.144, 0.346) 
8 = (0.740, 0.539, 0.330, 0.072, 0.104). 



The corresponding graph G is shown in Fig. 4. 

From the algorithm, we find the Hamiltonian circuit, P1P2P3P4P5, 
corresponding to the permutation tt = (12354). Therefore, we have the 
values shown in Table II. 

The routing strategy is to use route 777 for the fraction a, of the time. 

Example 2: There are also five paths between A and B. However, 
the desired traffic vector x' and the blocking probability vector p' are 
slightly different from those in Example 1. 

x ' = (0.191, 0.231, 0.220, 0.072, 0.242) 



Thus, 



p' = (0.7, 0.7, 0.6, 0.6, 0.3). 

5' = (0.638, 0.770, 0.550, 0.144, 0.346) 
8' = (0.445, 0.539, 0.330, 0.072, 0.104). 

The corresponding graph G' is shown in Fig. 5. 

From our algorithm, we find the Hamiltonian circuit, P 2 P\PzP5Pa, 
corresponding to the permutation it' = (21354). Therefore, we have 
the values shown in Table III. 

It is easily verified that 

r = I a/F(m'). 
»'=i 




Fig. 4 — Corresponding graph G for Example 1 . 
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Table II — Values of coefficients 
for Example 1 



Coefficient 



Route it, 



1 P,P 2 P 3 Pr,P4 0.897 

2 P-,P:iPt.P*Pi 0.032 

3 P;,Pr,P,P 1 P2 0.012 

4 P4P1P2P3P5 0.042 

5 PsP^PiPaPa 0.017 



Note that there is another Hamiltonian circuit in G', namely, 
PiP 2 PzPbPi, which gives an alternative cyclic routing realization as 
shown in Table IV. 

Again, it is easy to verify that 



a 



x' - £ afFGrf ). 



■1 



IV. CYCLIC APPROXIMATIONS 

As we mentioned earlier, the desired traffic vector x determined by 
the linear programming module can perhaps not be realized by a cyclic 
set of routes. In that case, we provide a routing strategy for approxi- 
mating x by modifying our cyclic routing algorithm. This is most easily 
explained in terms of an example (this one was taken from data 
generated by a 28-point simulation of Cardwell 6 ). 

In this example, there are 8 paths from A to B. Table V shows the 
appropriate data. 

Note that path 1 assumes the full traffic load, i.e., X\ = 1 - pi, which 
can be achieved if and only if every call requested is first attempted on 
path 1. In fact, this is a typical case of the existence of a least expensive 
direct line between two cities in a large toll switching network. 




Pa '3 

Fig. 5— Corresponding graph G' for Example 2. 
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Table III — Values of coefficients 
for Example 2 

Coefficient 
i Route ir't a,' 

1 P,PzP 5 P 4 P 2 O103 

2 P 2 P,P 3 Pr,P 4 0.730 

3 P 3 P 5 P,P,Px 0.109 

4 P4P2P1P3P5 0.042 

5 P5P4P2P1P3 0.017 



There are several reasons why this x cannot be realized by cyclic 
routing. For example, 

(i) Path 1 assumes the maximum possible traffic load, i.e., X\ = 
1 — pi, this cannot happen with cyclic routing. 

(ii) Traffic flow is highly unevenly distributed; in particular, paths 
3, 6, and 8 get no traffic at all. 

(iii) £, Xi 9* 1 - pip 2 •■- Pa. 

Let us form the graph G as described in the cyclic routing realization, 
namely, G has vertices {Pi, P 2 , • • • , Ps) and there is a directed edge 
from Pi to Pj if Si < oj (see Fig. 6). 

Here, G has no Hamiltonian circuit (and, in fact, is not even 
connected). In this case, we approximate x by taking a combination of 
(possibly trivial) disjoint circuits in G. The precise way this is done is 
described by the following algorithm. 

Algorithm for approximate cyclic routing 

(i) Calculate 

Xk . . PkXk , . 

Ok = and 8 k = , 1 < k < n. 

l-p k l-pk 

(ii) Relabel the a*, if necessary, so that a t > 02 > • • • > o n . If o„ = 
0, define t by a t > = o t +\, if a„ — 0. Otherwise, define t to be n. 

(iii) Set tti <— 1, 7 <— 1, i <— 2, ^ <— 5i, 2 <— 1. 

(ii>) If ^ > ff», go to (w). If v < a,, insert 1 after z in the cycle 
representation of vj. 



Table IV — Alternative values of 
coefficients for Example 2 



Coefficient 



Route 7t" 



1 PiPiPaPsP* 0.591 

2 P2P3P5P4P1 0.339 

3 PaPsftP.ft 0.012 

4 P 4 Piftp3Pa 0.042 

5 Pr,P*PiPzP a 0.017 
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Table v — Example not realized by cyclic routing 













Approxima- 




P 


X 


a 


8 


tion 


1 


0.30696 


0.69304 


1.00000 


0.30696 


0.69304 


2 


0.23850 


0.00084 


0.00111 


0.00026 


0.01832 


3 


0.30935 








0.00233 


4 


0.45325 


0.07555 


0.13818 


0.06263 


0.06989 


5 


0.28685 


0.13343 


0.18710 


0.05367 


0.12343 


6 


0.33891 








0.00116 


7 


0.60274 


0.09685 


0.24378 


0.14694 


0.08959 


8 


0.48781 








0.00196 



(v) liy> Si, set y <- 6\, z «- i. If y < 6\, y and z are unchanged. If 
i < t, go to (iv). If i = t, go to (uii). 

(vi) Set j *-j + 1, vj *-i,y*~ ^, z <- i. If i < t, set i *- i + 1 and 
go to (iv). If i = t, go to (mi). 

(uii) The routing strategy is given by using m, 7r 2 , • • • , ir a as follows. 

Let 





Pijt = Ok — 8irr>(*)j 




«,.* = =-j5- for £ £ 

Z. P'J 

J 


Use the route 





for the fraction 



of the time where 



Wi,i l W2,i 2 • • • fla,!. 



ai.i,a2.i, • • • «s,i s 



^A.i* = [**> W*(j*), 7T* (j*), • • •] 




Fig. 6— Corresponding graph G for the Cardwell example." 
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Continuing the example, we apply the above algorithm to the values 
in Table V. This results in the permutations -n\ = (1), tt 2 — (547), and 
7T 3 = (2). The corresponding as are given in Table VI. Note that we 
only use 5 of the 8 paths and the total blocking probability is 0.0057. 
If the blocking probability turns out to be too high to meet the required 
grade of service, we can make use of the remaining three paths in 
carrying the overflow by the modification shown in Table VII. The 
traffic flow generated by this routing strategy is listed in Table V 
under Approximation. Note that the approximation to the desired 
traffic flow is quite good. 

As we pointed out before, one of the main reasons we cannot achieve 
the desired traffic flow exactly is that this is inherently impossible to 
do using convex combinations of the available routes, because of 
premature termination or inadequate constraints in the linear program. 
In the next section, we examine a method for correcting this difficulty. 

Table VI — Values of coefficients for 
the Cardwell example 



Route Coefficient 

P,P 5 P<P7P2 «iia 2fl a 3 2 = 0.13132 

P,P 4 P7P5P2 a,.,a 2 .4a3.2 = 0.27634 

P,P 7 P 5 P«P2 a1.1aa.7a3i = 0.59234 



VII. UPPER BOUNDS 

Again, we consider a set of paths Pi, P2, • • • , Pn connecting two 
points and having blocking probabilities pi, P2, • • • , p n , respectively. 
The traffic flow on path P, cannot exceed its capacity, namely, 1 — />,. 
Thus, an immediate upper bound on Xi for any realizable traffic flow 
vector x is 

Xi ; < 1 — pi for all i. 

Similarly, for any two paths P, and P,, the total amount of traffic they 
can carry is 1 — piPj. Thus, if x is realizable then 

Xi + £/<!- PiPj. 



Table VII — Values of coefficients for 

modified routes of the Cardwell 

example 



Route Coefficient 



P,P 5 P 4 P7P2P6P8P3 0.13132 

P,P4P7P 5 P2P3P6P8 0.27634 

P,P7P5P4P2P 8 P3P6 0.59234 



1 784 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1 981 



More generally, for any set of k indices i h • • • , ik, if x is realizable then 
Xit + ••• +x ik <l-pi 1 ••• p* r (5a) 

Furthermore, for k = n, eq. (5) must hold with equality, i.e., 

Xi + • • • + x n » 1 - Pi • • • Pn. (5b) 

It is interesting to note that conditions (5a) and (5b) are also 
sufficient conditions for the realizability of x as a convex combination 
of flow vectors F(tt). The proof is not difficult. Basically, it is as follows. 
Suppose x is an extreme point of the polytope p defined by the 
intersection of the half planes (5a) and the hyperplane (5b). Then x 
must satisfy at least one of the equalities in eq. 5a with equality, say 
without loss of generality. 

Xi+ ••■ + x r = l-pi ••• Pn r<n. (6) 

We now use induction on n and express the point x' = (x\ • • • x r ) as 
a convex combination of the r! flow vectors associated with paths 
Pi, • • • , P r . We next consider all the inequalities in eq. (5) which 
contain *i, • • • , x r as well as other x's. Typically, we might have 

Xi + • • • + X r + X h + • • • + Xj m < 1 - Pi • • • PrPj, ' ' • Pu 

By eq. (6) 

X h + • ■ • + Xj a </?!••• Pr(l ~ Pi, • ' ' Pj)- 

Again, we can use induction, this time on the new variables yk = 
Xklp\ • • • Pr, r < k < n, which satisfy the required analogues of eqs. 
(5a) and (5b). Finally, we piece together these two convex combinations 
to get the desired representation for x. Since p is convex, then we are 
finished. 

Of course, in actual practice some appropriate subset of the ine- 
qualities in eq. (5a) would be used in the upper bounding process (see 
Ref. 1). 

VIII. CONCLUSIONS 

In this paper we give necessary and sufficient conditions for deter- 
mining whether a desired traffic flow vector (as specified by the linear 
programming solution portions of the algorithm) can be realized from 
a cyclic set of routes. The algorithm to verify the necessary and 
sufficient conditions can be implemented in 0(n log n) time. When the 
conditions are not met, we propose an approximation method which 
uses several smaller cycles rather than a single cyclic set of routes. 

In connection with the results we described earlier, it would be of 
interest to know what proportion of the volume of the polytope 
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spanned by the n! flow vectors F(rr) can in general be reached by cyclic 
routes. For n = 3, it seems that we can always cover at least % of the 
volume (actually, area in this case; see Figure 3). We have not yet 
examined the general case. In fact, we do not even know whether or 
not cyclic routes always span a positive fraction of the volume (inde- 
pendent of n). 
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